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Abstract 

We propose a simple computational procedure in order to resolve the degen¬ 
eracy, which invariably exists on the background of fluid motion associated 
with a channel of inhnite extent. The procedure is applied to elucidate the 
bifurcation structure for the particular case of laterally heated flow with the 
addition of a perturbative Poiseuille flow component. The introduction of a 
symmetry breaking perturbation as the simplest imperfection alters the bi¬ 
furcation tree of the original shear flow. As a result, the previously unknown 
higher order nonlinear solutions for the unperturbed flow are discovered, 
without implementing classical stability theory. 

Keywords: Bifurcation, nonlinear stability, 2:1 resonance, homotopy, 
symmetry breaking perturbation 
PACS: 47.10. -g 


1. Introduction and Methodology 

It is well known in quantum mechanics that a spectral line composed by 
different components can be split under the presence of a static perturbative 
magnetic held. This fact, the so-called Zeeman effect, implies that multiple 
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orbits of an electron can be generated with the introdnction of a pertnrbation. 
Similarly in this stndy, we show that degenerate states of shear flow can be 
detected with the introdnction of a compntational pertnrbation. 

Tnrbnlent channel flow is not only one of the fnndamental themes pnrsned 
in the held of hnid physics, bnt also of great engineering importance. Histor¬ 
ically, attempts to compntationally model tnrbnlent how and nnderstand its 
origins, strnctnre and evolntion were based on statistical approaches. Since 
the details of the aperiodicity (in velocity helds for example) observed in ex¬ 
periments or in nnmerical simulations are not of prime interest, one usually 
tries to characterise turbulent systems by their statistical and time averaged 
properties. Although it is widely accepted that the Navier-Stokes equations 
should provide the correct basis for an accurate description of turbulent hnid 
how, it is regrettable that rather little information from the equations, enters 
into the statistical analysis. 

More recently, the sequence of bifurcations/deterministic approach has 
emerged as a popular tool due to the discovery of coherent solutions that 
are universal. They exist in a variety of hnid hows as invariant sets within 
the turbulent regime and they have been obtained without the presumption 
of favourable or pre-set conditions deduced from empirical laws. Moreover 
the discovery of coherent structures indicates the existence of order within 
the evident aperiodicity of turbulence, that in turn could be used to probe 
aperiodic motion in a deterministic way. Therefore the sequential bifurcation 
approach, that identihes coherent structures and examines their stability, can 
potentially overcome many difficulties associated with transitions to turbu¬ 
lence. 

In order to describe in general terms the sequential bifurcation approach 
that will be followed in the next sections, we introduce the governing equation 
of hnid motion: 

-^ + L{R)u + N{u)=0. (1.1) 


In the set of partial diherential equations (1.1) u describes the deviation from 
the basic (undisturbed) state of the system, L{R) is a linear and N{-) is 
a nonlinear operator that involve partial derivatives. In (1.1) R represents 
collectively the parameters of the system, such as the Prandtl, the Rayleigh 
or the Reynolds numbers. In our sequential bifurcation approach we will 
be employing (1.1) and we hrst obtain the stability of the basic state via 
the introduction of inhnitesimal disturbances in (1.1) and by ignoring the 
nonlinear terms. The nonlinear states grow from the stability boundary of 
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the basic state, that is typically characterised by f{R) = 0, where / is the 
neutral surface of the basic state. If we repeat the aforementioned procedure 
of perturbation on the obtained nonlinear states, we could hnd a successive 
sequence of bifurcations for the higher order states en route to turbulence. 

We intend to deduce, or identify, numerically coherent structures for given 
values of the control parameters, for the case when the fluid is interacting 
with plates of inhnite extent. This is because the flow, in this case, is free from 
the influence of the presence of end boundaries at inlet or the outlet in the 
upstream and donwstream ends of the channel. Additionally the proposition 
of the inhnite channel allows the introduction of Fourier expansions into the 
simulations for the streamwise and spanwise directions. 

From the theoretical viewpoint, let us suppose a numerical channel with 
the gap width, 2h, under the periodic boundary conditions with wave num¬ 
bers a and (3 in the streamwise and spanwise directions.This channel is the 
subspace of another numerical channel with the same gap width but with 
(subharmonic) wave numbers and where n,m are arbitrary positive 
integers. An exact equilibrium state realised in the fundamental channel is 
dehnitely an exact state in the subharmonic channel.The latter exact state 
may bifurcate from the former exact state, which is traditionally termed as 
the subharmonic bifurcation. Taking an inhnite number of the combinations 
(m,n) into account, multiplicity of the exact states exists inevitably in the 
channel with inhnite extent. In this manuscript we will provide an example 
of how to probe previously unexplored subharmonic branches by a simple 
perturbation, without the traditional approach of linear stability analysis 
using Floquet parameters. 

In the following section, we formulate the problem and identify the basic 
state determined uniquely from the homogeneous external conditions. We 
obtain the stability boundary of the basic state in section 3, and proceed to 
solve (1.1), thus establishing the fully nonlinear solutions up to and including 
the tertiary level in the sequence of bifurcations approach to turbulence. We 
conclude this manuscript with a general discussion of our results and future 
work. 


2. Formulation 

We consider an incompressible Boussinesq Newtonian huid bounded be¬ 
tween two vertical parallel plates of inhnite extent with diherent temperatures 
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Figure 2.1: Geometrical configuration of the problem. The antisymmetric 
prohle for pure Laterally Heated Flow (LHF) deviates from the exactly an¬ 
tisymmetric form under a small perturbation induced by the mean pressure 
deviation induced in the streamwise direction by the Poiseuille flow compo¬ 
nent. 


Tq-I-AT and Tq—AT, subject to a pressure gradient deviation from the hydro¬ 
static pressure, see Fig. |2.1[ The origin of the Cartesian coordinates system 


is positioned in the midplane of the fluid layer of width 2d, taking x, y, z as 
the streamwise, spanwise and wall normal directions with unit vectors i, j, 
k. We assume here the state satishes the periodic boundary conditions with 
wave numbers a and f3 for the x and y directions, respectively. The velocity 
V and the temperature T are governed by the following equations: 


d 

—V+V■Vv 
ot 

1 9 

= — 'Vti + g'j T t + u'V^v, 

P 

(2.2) 

|-T + n-VT 
at 

= 

(2.3) 

V • V 

= 0, 

(2.4) 
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where we have incorporated the hydrostatic pressure term g'jTox in vr. Here p 
is the density, k is the thermal diffusivity, v is the kinetic viscosity, 7 is the co¬ 
efficient of thermal expansion and g is the acceleration due to gravity. For the 
non-dimensional description of the problem we use d, (P/v and AT/ PrGr, as 
the units of length, time and temperature, respectively. The physical proper¬ 
ties of the system are characterised by the three non-dimensional parameters 
Gr = gjATP/P, the Grashof number that gives the strength of the heat¬ 
ing, R = Umaxd/i', the Reynolds number that measures the strength of the 
applied pressure gradient in the streamwise direction {Umax is the speed of 
Poiseuille laminar flow at the origin) and the Prandtl number Pr = vj k. 
Taking the non-dimensional form of the velocity, temperature and pressure 
deviations from the basic state as w, Q and H, we obtain the non-dimensional 

—'VU + Oi + 'V'^u, (2.5) 

( 2 . 6 ) 

0. (2.7) 


form corresponding to (2.2)-(2.4), 


d 

—u -I- u ■ vu = 

ot 

d 

-e + u-we = 

ot 

V • u = 


Hereafter, we restrict the scope of our study to the case of Tr = 0 and two- 
dimensional flow (no spanwise dependence of the flow), so that our primary 
motivation in the present study, the retrieval of resonant modes, via the in¬ 
troduction of symmetry breaking perturbative flows, may be realised simply. 
At the vanishing limit of Pr, the conductivity of temperature is predominant 


against inertia, the energy equation (2.6) decouples from the momentum 


equation (2.5) and the temperature variation becomes simply proportional 


to 2; independently from u. The basic solution of equations (2.5)-(2.7), which 


satishes the no-slip condition for the velocity {u = 0) on the boundaries at 
2; = ±1, is given by: 


Ub(:) = U„(z) i, U„(z) = (Gr/6)(z^ - z) + R(1 - 


( 2 . 8 ) 


For the conhguration of Figure 2.1 Squire’s theorem is applicable (see Squire 


( 1933[ )). We further assume that the pressure gradient remains in the mean 
unchanged by the amplitude of the secondary flow and we show that there 
exists degeneracy that is invariably linked with the possibility of inhnite 
streamwise length. In order to obtain the fully nonlinear equilibrium solution 
of these two-dimensional rolls we introduce the stream function p, a function 
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of (x, 2 :)-variables, with the velocity fluctuations u satisfying: 

= V x Ip 3 - 

The boundary conditions for ip are given by: 


/ n 

iP= — = 0 

oz 


= 0 at z = ± 1 . 


(2.9) 


( 2 . 10 ) 


Two-dimensional steady equilibrium roll solutions were therefore sought nu¬ 
merically employing the Chebyshev Gauss-Lobatto collocation point method 
combined with the Newton-Raphson iterative method via the following har¬ 
monic expansion for ip: 


N M 

1p{x,Z,t) = '^ ^ ttnm 

n=0 m=—M, nif^O 


(i-.r 


r„(.-) 


^im a(x—ct) 


( 2 . 11 ) 


where N and M are the truncation levels for the complex coefficients anm, 
and Tn{z) is the u-th order Chebyschev polynomial. We have incorporated 
the phase velocity c in the expansion, so that we can examine additionally 
the case, where a perturbative Poiseuille flow component is allowed to be 
present in the configuration of Figure |tH via the implementation of i? 7 ^ 0 
in (2.8). In this case calculations are performed on a moving frame that is 


phase locked with the nonlinear solutions bifurcating from the neutral curve 
boundaries. The factor (1 — has been inserted in the expression for ip in 


order to take into account the boundary conditions expressed by eqs. ( 2 . 10 ). 

By subs titut ing the expansion (2.11) into the streamwise projection of 
the curl of (2.5), multiplying the resulting equations by 


dxe 


llOLX 


I = 


1 , 2 ,..., we can evaluate the expansions ( 2 . 11 ) at each collocation point, i.e 


we obtain a set of nonlinear algebraic equations for the complex coefficients 
anm- Further details of the numerical method employed here to obtain nonlin¬ 


ear solutions have been presented recently in Generalis and Fujimura (2009) 


and Itano and Generalis (2009); Generalis and Itano (2010) in relation to 


different types of shear ffows. We elaborate on the results of our simulations 
in the next section. 


3. Results 

3.1. Linear stability 

We dehne the neutral stability curve, as the curve that provides the 
critical Gr number as a function of the wavenumber a at a pre-specihed 
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value of Reynolds number say, Grc{a,R). In Figure 3.2, the neutral curve 
Gr = Grc{a, R = 0), is given as a thick solid curve on the a-Gr plane. The 
thin solid curves on the left hand side of the a-Gr plane are also neutral 
curves, but they depict the Grc{na,0) relationship in the case, where the 
scale of a has been reduced by a factor 1/n, n = 2, 3,4. These are described 
by the functional relationship Gr = Grc{na) and we term them subharmonic 
neutral stability curves. At a given Gr> Grc{na,0), an inhnitesimal distur¬ 
bance added onto the basic flow may grow exponentially. The eigen function 
generated from Grc{a, 0) is referred to as the fundamental mode, while those 
from the duplications Grc{na, 0) {n = 2, 3,4) are referred to as subharmonic 
modes. 


Table 3.1: Wavenumber and Grashof number at bi-critical point, (ccc.n, Grc,n), 
at R = 0 are listed, which are determined as the solution of Grc{a,0) = 
Grc{na,0). At the bi-critical point, the nonlinear branch may be generated 
from a certain noninear interaction between the fundamental and subhar¬ 
monic modes. For example, the branch T shown in Fig. |3.2| originates at a 


second bi-critical point, (q:c, 2 , listed in Tab,3.1 


n 


Grc,n 

2 

0.8631 

604.3 

3 

0.6166 

791.9 

4 

0.4761 

999.9 

5 

0.3879 

1212.7 

6 

0.3277 

1426.1 

7 

0.2839 

1639.1 

8 

0.2506 

1851.9 


3.2. Resonance in the infinite extent case 

Since we assumed that the system has an inhnite extent, it is necessary to 
keep in mind the fact that the characteristic wavelength is selected by nature, 
that is, the wavelength dominating the bifurcating nonlinear state should not 
be pre-assumed by our calculations. It is possible that any combination that 
consists of a number of eigenmodes, associated with the neutral curves, can 
develop into some resonant nonlinear state in the system with inhnite extent. 
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Figure 3.2: The neutral curves for the fundamental (solid) and subharmonic 
modes (dashed) at i? = 0, which is described as Gr = Grdna, 0) n = 1, 2, 3,4. 
The critical Grashof number is Gr = 495.6 at the wavenumber a = 1.344. 
The intersections of the curves, which are termed as “bi-critical points”, are 
all candidates of the origins for the nonlinear states generated by a certain 
nonlinear interaction between the different modes, the pure (nonlinear states 
from the neutral curve) and the subharmonic neutral curves. This is the so 
called “resonance”. 


The case of the conventional Rayleigh Benard convection is an example, 
where, square or hexagonal cells are constituted by the combination of two or 
three self resonating eigenmodes. In Table 3T, we list the specific parameter 
values {a, Gr), for the “bi-critical points”, at which the neutral stability curve 
intersects the subharmonic stability curves, Grc{a) = Grcijia) {n > 1). At 
these parameter values the fundamental mode may interact the subharmonic 
modes, so that a non-linear state can be established. 









a 


Figure 3.3: Nonlinear branches at {Gr, R) = (1100,0). The thick curves 
marked by F± is the secondary branches bifurcated directly from the basic 
state on the neutral curve, Gr = Grc{a,0). In particular, the curve marked 
F_ connects to those of F 4_3 and F 3_2 sequentially on the thin curves with 
S„, which refer to the n-th subharmonic duplications of the F+. The gray 
curve with symbol T and the other gray curves depicted in the left hand 
side of the graph correspond to the travelling wave solutions obtained in the 
present procedure. The lower case corresponds the location in the graph 
where the corresponding flow pattern illustrated in the Fig]3.4|are obtained. 


3.3. Previous work 

Previous studies on the subject have remained an ambiguous point in 
their attempts to clarify the nature of the interaction that the fundamental 
and subharmonic modes play, as well as the crucial role that these inter¬ 
actions play for a system with inhnite extent as a whole. For the case of 
i? = 0, as in the present model, Nagata and Busse (1983) calculated equilib¬ 
rium states, which should be expected to bifurcate from Grc{a,0), but the 
hnite-amplitude equilibrium states were obtained only for a narrower range 
than expected. Let us explain this. Suppose that the laminar state looses 
its stability at a given Gr against a perturbation with the wave number a 
satisfying a_ < a < a+, where a± are the solutions of Gr = G'rc(a,0). 
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Figure 3.4: Contour plots for stream function for Gr = 1100, R = 0. The 
contour levels are indicated in light grey scale for small and in dark grey 
scale for large values and for: (a) a = 1.60 on F+, (b) a = 0.98 on F+, (c) 
a = 0.42 on F_, (d) a = 0.48 on F_, (e) a = 0.64 on F 4 _ 3 , (f) a = 0.88 on 
F 3 _ 2 , (g) a = 0.42 on T, (h) a = 0.88 on T. In each plot the state between 
the boundaries is depicted for one wavelength in the ordinate. The location 
(a)-(h) of the corresponding flow pattern in the {a, Gr) plane is indicated in 


Fig,3.3 


Fujimura and Mizushima (1987) claimed that the equilibrium state can exist 


only in the range of wave number a+/2 < a < a+. Fujimura and Mizushima 


(1987) could partly explain the mystery mentioned above. Employing weakly 
nonlinear analysis, they showed the possibility that an interaction of the fun¬ 
damental mode with the subharmonic mode may constitute a nonlinear equi¬ 
librium state with subharmonic wave number, which could not be obtained in 


Nagata and Busse (1983). The case of Rayleigh-Benard convection was also 


reported in Busse and Clever (1974) and Mizushima and Fujimura (1992). In 


Mizushima and Fujimura (1992) the discrepancy in the bifurcation diagram 


of Nagata and Busse (1983) was accounted for, by attributing one of the 
’hidden’ states as the state generated by the mutual interaction between the 
fundamental and subharmonic modes n (= 2, 3,4, • ■ ■). This was termed as 
1 : n resonance, according to them. 


3.4- Unravelling the mystery of the resonant region 

By introducing to the system with inhnite extent a perturbation breaking 
the symmetry of the basic state, we will obtain possible resonant nonlinear 
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states at Gr = 1100 and R = 0 (see Fig. 3.3 (a)). Following Herbert 


(1980), we characterise states in terms of the kinetic energy of the flnctnation 


normalised with respect to the basic flow, E, 


m=oo 

e=Y,e„ 

m=l 


Em = ^ 


r- + l 


1-1 


^‘I'K j Ot 


N 


IVV’mP dx dz, = + c. c. 


n=0 


(3.12) 

In Fig. 3.3, the secondary branches bifnrcated from the nentral cnrve Grc{a, 0) 
at a = Q!± for Gr = 1100 are depicted as thick cnrves with symbol F±. F± 
stands for the secondary state that bifnrcates from the fundamental nentral 
cnrve at a±. The three thin cnrves, S„, are dnplications with scale factors n 
(= 2,3,4) from F+. S„ stands for the secondary state that bifnrcates from 
the subharmonic nentral cnrve at a±/n. T stands for a tertiary (travelling 
wave) state bifnrcated from the secondary states. 

Branch F+ bifnrcated from a+(= 1.922) terminates at the node on S 2 , 
a = 0.973. Along the F+ as a decreases, the flow pattern changes from single 
vortex state (Fig. |3.4K a)) to two vortex state (Fig. |3.4K b)). On the other 
hand F_ bifnrcated from a_ (= 0.430) terminates at a = 0.481 on S 4 . Along 
the F_ as a increases, the flow pattern changes from single vortex state (Fig. 
|3.4K c)) to fonr vortex state (Fig. 3.4[ d)). The branch F_ joins with a branch 
F 34 at a = 0.481 on S 4 , and F 34 connects to a branch F 23 at a = 0.644 on S 3 . 
Fmn; rn,n integers, refers to the nonlinear state bifnrcated from the nentral 
cnrve that exists in the m, n resonant region. At the node on S„ n vortex 


state (Fig. 3.4'd)-(f)) is observed. Finally, F 3_2 terminates at the node at 
S 2 . There is only the branch S 2 between a = 0.888 and 0.976. Note that the 
first Fonrier mode is absent for S 2 . This gap is the mysterions region where 


no state could be found in the previous study Nagata and Busse (1983). 


3.5. Introduction of perturbation i? 7 ^ 0 

In the following, we will show a hidden non-linear branch, both ends of 
which connect to the non-linear branches emerging from the neutral stability 
curve G'rc(a, 0) and the duplications Grc{na, 0) (n = 2, 3,4). Here, we intro¬ 
duce a perturbation, a small plane Poiseuille flow component, to the system 
with i? = 0, by employing a non-zero value of R. As the degenerate spectra 
are uncovered by the Zeeman effect, it is expected that the introduction of 
such a perturbation consisting of a net drifting in the streamwise direction 
will induce a structural instability in the system with i? = 0 , which origi¬ 
nally satisfies a relatively higher symmetry. Although the introduced value. 
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i? = 1 , seems to be subtle, the resulting bifurcation diagram is substantially 
richer than that obtained by means of the classical step-by-step procedure 
for the sequential bifurcation. It is of particular interest that the uncovered 
nonlinear state in the diagram at i? = 1 survives even when R = 0. 


3.6. The discovery of the branch T 

For example, branch T in Fig j3.3| corresponds to a nonlinear state uncov¬ 
ered in this procedure. The branch T ’hlls’ the aforementioned gap, which 
has been previously unexplained properly. This was found systematically 
by the present procedure utilising the structural instability of the system 
with i? 7 ^ 0. We should emphasize that, even at i? = 0, the branch T pos¬ 
sesses a non-zero phase velocity, c, which originates from the system that 
was perturbed with R = 1. Instead if we select the opposite sign of the 
perturbation, R = —1, in the procedure, the same kind of solution drifting 
in the opposite direction can be obtained. Therefore two drifting solutions 
have been identihed via homotopy followed by its reversal. Extending over 
the previously unexplained gap, the branch T forms a bridge between the 
branches, F+ and 83 - 2 . The contour plot of T is given in Fig j3.4[ g) and (h), 
which are somewhat reminiscent of the vortices in Fig j3.4K c) and (f). This is 
probably relevant to the existence of the mixed mode suggested by Fujimura 
and Mizushima Fujimura and Mizushima ( 1987[ ), even though they obtained 
evidence of its existence only at a few sets of parameters by weakly nonlinear 
analysis where the obtained results would be only valid in the vicity of the 
neutral curves. 


3. 7. Remainder of nonlinear branches 

Our numerical simulations have revealed several more nonlinear branches, 
which are shown on the left hand side of Fig j3.3[ At hrst, there are additional 
branches in the region 0.29 < a < 0.59, see Fig j3.3[ We have only depicted a 
sample of tertiary states to keep the hgure uncluttered. The branch S„ serves 
as the platform for bifurcation of three higher order distinct branches in the 
regions, 0.467 < a < 0.510, 0.316 < a < 0.418, and 0.318 < a < 0.333. 
Moreover, another tertiary state was found in a tiny region a ~ 0.4712 — 
0.4873. All the branches indicated by thick curve including T possess the 
phase velocity, which are described in the Figj3.5[ 
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Figure 3.5: Phase velocity for Gr= 1100 and R = 0. 


4. Concluding remarks 

4.I. Numerical validation of obtained result 


In Table 4.2 we show a comparison between our nonlinear results and 
those of Mizushima and Saito (1987) for the value of Ei of (3.12). Despite 
the different truncation levels employed in the two studies, the results of the 
corresponding simulations agree to a high degree of accuracy. 


4-2. Subcriticality 

Additionally, we would like to note a subtle change in bifurcation diagram 
on a — Gr plane occuring at the wavenumber a+ between Gr = 900 and 
1100. The bifurcation of the nonlinear branch, F+ from the basic state is 
supercritical at Gr = 900. However, beyond a value of Gr between 900 < 
Gr < 1100, the bifurcation of nonlinear state F_|_ becomes subcritical. This 
would be better shown in Fig j3.3[ where the bifurcation is subcritical, as the 
existence of region of the nonlinear state moves towards larger wavenumber 
in the hrst instance. 
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Table 4.2: Comparison of nonlinear results with Mizushima and Saito (1987) 
for the energy Ei of (3.12) of the pure mode for various {Gr,a) values. In 
the present work the truncation level (L, M) = (21,18) has been used for 
all nonlinear calculations. In Mizushima and Saito ( 1987| for the nonlinear 
calculations L = 40 was employed with the value of M as indicated. 


Gr 

a 

Mizushima &: Saito (M) 

Present results (M = 18) 

500 

1.25 

2.41706 X 10-^ (4) 

2.41746 X 10-^ 

500 

1.30 

2.53302 X 10-3 (4) 

2.53306 X 10-3 

600 

0.90 

9.25916 X 10-^ (6) 

9.25916 X 10-^ 

600 

1.30 

6.14480 X 10-2 (6) 

6.18481 X 10-2 

2000 

2.00 

1.88510 X 10-1 (12) 

1.88509 X 10-1 

5000 

1.30 

3.70190 X 10-2 (16) 

3.70657 X 10-2 


4-3. Topological change of bifurcation diaggram with increase of Gr 

We explain now the topological change of the bifurcation diagram with 
decreasing of Gr. Although it is not illustrated here, for Gr = 500 around 
the critical Grashof number, the bifurcation of the nonlinear branch is su¬ 
percritical at a± and has the maximum value in at a ~ (a+ -|- a-)/2, as 
it is a convex one single curve. The nonlinear branch does not have phase 
velocity (stationary nonlinear solution). 

This is not the case for higher values of Gr. In our study we have found 
that the single curve of the nonlinear state ’breaks-up’ over the bi-critical 
point, that is, Gr > Grc^ 2 i due to the existence of ’’resonance”. At the 
bi-critical point where the disconnection occurs, the two split branches are 
connected with the hidden resonant neutral curve, T. Note that the tertiary 
branch T exists at smaller Grashof numbers than Greg around S 2 . 


f-f- Topological change of the bifurcation diagram from R = 1 to R = 0 
Initially the pure mode nonlinear solution, which bifurcates from the 
neutral curve at a = 1.9223 (see Figure 3.3) has a phase velocity (c = 


0.61447, a = 1.9223), (see also Figure 3.3 for comparison), due to the pres¬ 
ence of R, follows the path of the pure LHF, but now at a = 1.06961 the 
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pure mode (that bifurcates from the left hand side of the neutral curve) 
folds and continues within the resonant region, where it connects with the 
subharmonic 1 : 2 resonant mode at a = 0.8850. Similarly the solution 
that bifurcates from the left hand side of the neutral curve at a = 0.4299 
(with phase velocity (c = 0.48633, a = 0.42994), also connects with the 
1 : 2 resonant mode branch at a ~ 0.97251. Therefore the entire branch 
of the secondary state extends naturally over to the resonant region via two 
oppositely drifting waves with opposite phase velocities, with both drifting 
waves connecting with the 1:2 subharmonic (resonant) nonlinear solution. 


As shown in Figure 3.3, the phase velocities of these two secondary branches 
(including their loops) are ’bridges’ between the disconnected parts of the Ei 
curve and the 1 : 2 subharmonic (resonant) region for [Gr = 1100, R = 1) 
over a range of the wavenumber a. 

4-5. Advantage in this procedure : summary 

The measure of the influence of the PFC on the bifurcation spectrum of 
the laterally heated flow without PFC is shown in Fig. 3.3[a), if subsequently 
we gradually remove the imposed constant pressure gradient. The important 
aspect here is to note that our calculations have shown that the connecting 
secondary nonlinear states, extensions of the Ei curve for the case R = 1 
and embedded in the resonant region, remain if via homotopy we remove 
the added value of R and slowly arrive ’back’ to the case when there is no 
PFC, i.e. pure LHF. In other words, while for R = 1 there is topologically 
two 77—curves, that naturally connect the pure mode (bifurcating from the 
right and left hand side of the neutral curve) disconnected branches and the 
resonant bifurcation branches via the introduction of the oppositely moving 


drifting waves, as shown in Figure 3^, the topology of the bifurcation tree 
changes substantially for R = 0. The difference (when the constant pressure 
gradient has been removed, i? = 0), is that the connecting branches break¬ 
up now and the discontinuity points become bifurcation points of higher 
order, i.e. tertiary, states. The case for R = 0, contains a much richer 
structure than the case R = 1, since for the latter the higher order part of 
the bifurcation is completely absent and only the secondary nonlinear states 
are present. Without the introduction of the Poiseuille flow perturbation 
though, we would have not been able to discover the rich structure hidden 
in the case of i? = 0. 
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